Physisorption of an electron in deep surface potentials off a dielectric surface 
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We study phonon-mediated adsorption and desorption of an electron at dielectric surfaces with 
deep polarization-induced surface potentials where multi-phonon transitions are responsible for elec- 
tron energy relaxation. Focusing on multi-phonon processes due to the nonlinearity of the coupling 
between the external electron and the acoustic bulk phonon triggering the transitions between sur- 
face states, we calculate electron desorption times for graphite, MgO, CaO, AI2O3, and Si02 and 
electron sticking coefficients for AI2O3, CaO, and Si02. To reveal the kinetic stages of electron 
physisorption, we moreover study the time evolution of the image state occupancy and the energy- 
resolved desorption flux. Depending on the potential depth and the surface temperature we identify 
two generic scenarios: (i) adsorption via trapping in shallow image states followed by relaxation to 
the lowest image state and desorption from that state via a cascade through the second strongly 
bound image state in not too deep potentials and (ii) adsorption via trapping in shallow image 
states but followed by a relaxation bottleneck retarding the transition to the lowest image state and 
desorption from that state via a one step process to the continuum in deep potentials. 
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I. INTRODUCTION 

Image states, arising from the polarization-induced in- 
teraction between an electron and a surface, offer the 
possibility for electron trapping at a surface. Since 
their original prediction^ for the surface of liquid and 
solid He, they have been extensively studied for metal- 
lic surfaces.-^^ But image states also exist for dielectric 
surfaces provided the electron affinity of the dielectric is 
negative, that is, the vacuum level falls inside the gap 
between the valence and the conduction band. Image 
states are then the lowest unoccupied states and should 
hence allow for temporary trapping of external electrons. 
So far image states at a dielectric surface have been only 
observed for graphite but they are expected for other 
dielectrics with negative electron affinity as well, for in- 
stance, boron nitride^ and the alkaline earth oxides.^ 

Based on the idea of a two-dimensional electron sur- 
face plasmajiSr— electron trapping in image states has 
been suspected for a long time to be responsible for the 
build-up of surface charges at plasma walls. We have 
recently proposed therefore to consider the charging of 
a plasma wall as an electron physisorption procesSi ^^d^ 
Indeed, for plasma walls with negative electron affinity 
image states should contribute to the very beginning of 
the charging process, when the wall carries no charges 
yet and the image states thus fall inside the energy gap 
of the wall. Only with increasing surface charge image 
states are expected to play a less important role because 
the Coulomb barrier due to the electrons already resid- 
ing on the wall shifts image states to an energy range 
where they are destabilized by unoccupied bulk states. 
The later stages of charge collection most probably oc- 
curs via surface resonances or empty volume states.^® 

Regardless of its importance for charge collection at 
dielectric plasma walls, the electron kinetics in the im- 
age states of a dielectric surface is an interesting phe- 
nomenon in its own right. In addition, it is relevant in 
other physical contexts as well. For instance, (i) in elec- 



tron emitters, such as cesium-doped silicon oxide films 
with negative electron affinity, electron emission via im- 
age states reduces the operational voltage considerably^ 
(ii) in gallium arsenide based heterostructures surface 
charging can be used for the contactless gating of field 
devices,— and (iii) for the alkaline earth oxides, studied 
in the field of heterogeneous catalysisjiS"— the electronic 
surface states provide the environment for catalytic reac- 
tions. Some situations encompass electronic transitions 
from bulk to surface states, as it is the case for electron 
emitters, while for others, the electron does not penetrate 
into the bulk and the electron kinetics takes only place in 
surface states. Interesting questions in this case are the 
probability for temporary trapping in these states, the 
mechanism of electron energy relaxation at the surface, 
and the time after which a trapped electron is released. 

This is the concluding paper out of a series of three 
on the phonon-mediated physisorption of an electron in 
the image states of a dielectric surface. As in our pre- 
vious work, Refs. [2^ and [H (thereafter referred to as I 
and II), we investigate adsorption and desorption of an 
electron at finite temperatures assuming an acoustic lon- 
gitudinal bulk phonon to control energy relaxation at the 
surface. For the dielectric material we are considering, 
the level spacing of the lowest two bound states typically 
exceeds the Debye energy, implying that multi-phonon 
processes have to be taken into account. In I and II, we 
have studied desorption and sticking using an expansion 
of the energy dependent T matrix^"— allowing to calcu- 
late one- and two-phonon transition probabilities. This 
approach is however limited to very few materials, for in- 
stance, graphite and MgO. In the following we will adopt 
a different strategy, calculating multi-phonon transition 
probabilities due to the nonlinear electron-phonon cou- 
pling non-perturbatively. This allows us to calculate the 
desorption time and the sticking coefficient for the deeper 
surface potentials of CaO, AI2O3 and Si02. 

The remaining paper is structured as follows. In 
Sec. |TT] we briefly recall the quantum-kinetic approach 
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to physisorption. In Sec. IIIIl we calculate the multi- 
phonon state-to-state transition probabilities. In Sec. IIVI 
we present our results for the desorption time and the 
prompt and kinetic energy-resolved and energy-averaged 
sticking coefficient. In this section we also discuss the 
time evolution of the bound state occupancy and the 
energy-resolved desorption flux. Section |V] is devoted to 
the analytic treatment of a simplified two-state model, 
used to identify two generic physisorption scenarios into 
which we can classify the results of this paper as well as 
our previous results, before we conclude in Sec. I VII 



II. ELECTRON KINETICS 

As in and Il2i we describe the time evolution of the 
occupancy of the bound surface states with a quantum- 
kinetic rate equationi^i^ It captures all three character- 
istic stages of physisorptioni^i^ initial trapping, subse- 
quent relaxation and desorption. 

The time dependence of the occupancies of the bound 
states is given b y^^'^^ 



If the modulus of one eigenvalue, Aq, is considerably 
smaller than the moduli of the other eigenvalues, A^, 
a unique desorption time and a unique sticking coeffi- 
cient can be identified. In this case Aq governs the long 
time behavior of the equilibrium occupation of the bound 
states, n^'i ~ g-E^/ksTs ^ ^^^^j j^-g jj^yerse can be identified 

with the desorption time, X^^ = t^- In this case the 
bound state occupancy ri„(t) splits into a slowly varying 
part (t) given by the k = summand in Eq. Q and a 
quickly varying part n^n{t) given by the sum over k ^ 
in Eq. 

The "adsorbate", i.e. the fraction of the trapped elec- 
tron remaining in the surface states for times on the or- 
der of the desorption time, is given by the slowly varying 
part only, n^{t) — X]n"-ri(0- Differentiating nP^t) with 
respect to the time, 



^kinetic 



jfc(t)-Aon°(t) , 



(5) 



we can, following Brenig^, identify the kinetic energy- 
resolved sticking coefficient 



kinetic 



(O)p(O) 



nk 



(6) 



WknTlnit) + V TtWnkJk{t) (1) 



fe 



k 



= J2 Tnn'^W (t) + J2 nWnkJkit) , (2) 
n' k 

where Wn'n is the probability per unit time for a transi- 
tion from a bound state n to another bound state n' , Wkn 
and Wnk are the probabilities per unit time for a transi- 
tion from the bound state n to the continuum state k and 
vice versa and tj = 2L/vz is the transit time through the 
surface potential of width L, which, in the limit L — )> oo, 
can be absorbed into the transition probability. The ma- 
trix Tnm is defined implicitly by the above equation. The 
last term in Eqs. ^ and respectively, gives the in- 
crease in the bound state occupancy due to trapping of 
an electron in bound surface states. 

The probability for an approaching electron in the con- 
tinuum state k to make a transition to any of the bound 
states is given by the prompt energy-resolved sticking 
coefficient. 



prompt \ ^ TjT 



(3) 



Treating the incident electron flux as an externally 
specifled parameter, the solution of Eq. ([T]) describes 
the subsequent relaxation and desorption. It is given by 

Unit) = J2 f dt'e^-^'ei'') ^'i'^^^tWikjkit') , 

(4) 

where ei"'' and el'^'' are the right and left eigenvectors to 
the eigenvalue — A^ of the matrix T. 



giving the probability for both, initial trapping and sub- 
sequent relaxation. 

If the incident unit electron flux corresponds to an 
electron with Boltzmann distributed kinetic energies, the 
prompt or kinetic energy-averaged sticking coefficient is 
given by 



Ek^e:kke-^^''' 



(7) 



where /S^^ — kBT^ is the mean electron energy. 

The desorption flux, that is, the flux due to an electron 
that is not instantly reflected at the boundary but sticks 
to the surface for a flnite time can also be calculated from 
the occupancy of the bound surface states. FromEq. ([!]), 
we infer that the losses of the bound state occupancy 
increase the continuum state occupancy by 



"dT 



(8) 



As the electron remains in the surface potential for the 
time it needs to travel trough the surface potential the 
occupancy of the continuum state k is given by Uk = 
Ttrik- To obtain the energy-resolved desorption flux we 
multiply the occupancy of the continuum state k with 
the flux associated with the box- normalized state 

\(j)k)-'^'^ Thus, the energy-resolved desorption flux is given 
by 



]k{t) 



(9) 



which is well deflned in the limit L ^ oo. 
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III. TRANSITION PROBABILITIES 

The kinetic equations presented in the last section rely 
on the knowledge of the transition probabilities. They 
have to be calculated from a microscopic model for the 
electron-surface interaction. 

For a dielectric surface, the transitions are driven by 
phonons, whose maximum energy is, within the Debye 
model, the Debye energy hojo- Measuring energies in 
units of the Debye energy, important dimensionless pa- 
rameters characterizing the potential depth are 



En 

hojT 



and Ann' = 



En — En' 
huJri 



(10) 



where En < is the energy of the n*'' bound state. 

In I, we introduced the following classification for the 
potential depth. If —n -I- 1 > A12 > —n, we call the po- 
tential n-phonon deep. For the calculations in I and II, 
we considered only one- or two-phonon deep potentials, 
for which one- and two-phonon transition probabilities 
are sufficient. Dielectrics with two-phonon deep poten- 
tials, such as graphite or MgO, are however an exception. 
Many dielectrics, for instance, AI2O3, CaO, GaAs, or 
Si02 have more than two-phonon deep potentials. Hence, 
the more relevant situation is physisorption in deep sur- 
face potentials for which multi-phonon transition proba- 
bilities are required. 

To calculate multi-phonon transition probabilities for 
the one-dimensional microscopic model used in I and II, 
we briefly recall its main features. In short, for a di- 
electric surface, the main source of the attractive static 
electron-surface potential is the coupling of the electron 
to a dipole-active surface phonon*^ Far from the surface 
the surface potential arising from this coupling merges 
with the classical image potential and thus ^ 1/z. Close 
to the surface, however, the surface potential is strongly 
modified by the recoil energy resulting from the momen- 
tum transfer parallel to the surface when the electron 
absorbs or emits a surface phonon. Taking this effect 
into account leads to a recoil-corrected image potential 
~ l/{z + Zc) with Zc a cut-off parameter defined in I. 

Transitions between the eigenstates of the recoil- 
corrected image potential are due to dynamic perturba- 
tions of the surface potential. The surface potential is 
very steep near the surface. A particularly strong per- 
turbation arises therefore from the longitudinal-acoustic 
bulk phonon perpendicular to the surface which causes 
the surface plane to oscillate. 

The Hamiltonian from which we calculate the transi- 
tion probabilities was introduced in I where all quanti- 
ties entering the Hamiltonian are explicitly defined. It is 
given by 



ph 



H. 



dyn 

-ph 



where 



H: 



static 



(11) 



(12) 



describes the electron in the recoil-corrected image po- 
tential. 



(13) 



describes the free dynamics of the bulk longitudinal 
acoustic phonon responsible for transitions between sur- 
face states, and 



H 



dyn 
-ph 



E 



[q'\Vp{u,z)\q)cl,Cq 



(14) 



denotes the dynamic coupling of the electron to the bulk 
phonon. 

The perturbation z) can be identified as the dif- 
ference between the displaced surface potential and the 
static surface potential. It reads, after the transforma- 
tion Z ~> Z — Zr, 



Vp{u,z) = - 



z + u 



(15) 



where Ag = (e^ — l)/4(es -I- 1) with Cg the static dielec- 
tric constant. In general, multi-phonon processes can 
arise both from the nonlinearity of the electron-phonon 
coupling i?gZp;i as well as from the successive actions of 
^e~r>h encoded in the T matrix equation. 



rp _ frdyn 
^ ^e-ph 



H 



dyn 
e~ph 



GnT 



where Gq is given by 

Go^{E- Hf 



ph 



(16) 



(17) 



The transition probability per unit time from an elec- 
tronic state q to an electronic state q' encompassing both 
types of processes is given by^'^ 



27r 



\{s',q'\T\s,q)f 



x5{Es-Es, +Eq-Eq,) , 



(18) 



where j3s — {kBTs)~^, with Tg the surface temperature 
and \s) and |s') the initial and final phonon states. We 
are only interested in the transitions between electronic 
states. It is thus natural to average in Eq. (|T8)) over 
all phonon states. The delta function guarantees energy 
conservation. 

In our previous work, we have used an expansion of 
the T matrix to calculate multi-phonon transition rates. 
In principle this ensures that both linear and nonlinear 
terms in the interaction as well as successive actions of 
the perturbation are taken into account up to a certain 
order of the phonon process. However even for a two- 
phonon deep potential, taking all two-phonon processes 
into account is nearly impossible. The calculation be- 
comes feasible if two-phonon processes are only taken 
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into account for transitions not already enabled by a one- 
phonon process. This amounts to computing only the 
lowest required phonon order for a given transition, ne- 
glecting higher order corrections to it. For higher order 
phonon processes even this simplified strategy becomes 
unfeasible. A different approach is thus needed. 

From I and II we qualitatively know the relevance of 
the different types of multi-phonon processes for partic- 
ular electronic transitions. For continuum-bound state 
transitions, for instance, one-phonon processes are suffi- 
cient at low electron energies. We will therefore compute 
the transition probability between bound and continuum 
states in the one-phonon approximation. For transitions 
between bound states, we found that multi-phonon pro- 
cesses due to the nonlinearity of the electron-phonon cou- 
pling tend to be more important than the multi-phonon 
processes due to the iteration of the T matrix, unlike to 
what we found for bound state - continuum transitions 
(see I) or to what Siber and Gumhalter'^'^"'^^ found in 
the context of atom-surface scattering. Indeed, multi- 
phonon processes from the iteration of the T matrix give 
a minor contribution, unless resonances arising from the 
T matrix become relevant. This happens whenever the 
energy difference between two bound states is a multi- 
ple of the Debye energy. Resonances smoothen then the 
abrupt steps in the transition probability at the depth 
thresholds. Since the electronic matrix element between 
the first and the second bound state is the largest one, 
this effect is most pronounced for IA12I = n. 

In view of the above discussion we expect an approxi- 
mation which takes only the nonlinearity of the electron- 
phonon interaction non-perturbatively into account to 
give an acceptable first estimate for the multi-phonon 
transition rates. We denote this approximation the non- 
linear multi-phonon approximation. In particular, it 



should be sufficient for the identification of the generic 
behavior of multi-phonon-mediated adsorption and des- 
orption. 

Calculating multi-phonon processes due to nonlinear 
terms in the interaction potcntiaP^ amounts to a dis- 
torted wave Born approximation with the full interaction 
potential. Thus, the transition probability per unit time 
is given by 



e 



X 5{E,+Eq~Eq, -E,,) 



(19) 



To evaluate the multi-phonon transition probability, 
we use H'^I.pfi in the form of Eq. (IT5|) . The transition 
matrix element in Eq. (I19p is then given by 

dz0;(z) [v(z + u) - v(z)] cl>g,{z)\s') , (20) 



where v{z) = — (e^Ao)/^. Introducing dimensionless 
variables x ~ z/aB, the Fourier transform of the static 
potential 



(21) 



(22) 



and the state- to-state matrix element 

J Xc 

the transition probability can be rewritten as 



2tt 



E 



-PsE^. 



dp f°° dp 



2tt 



2tt 



Vip)v*ip)fqq'ip)fnn'{p){s 



W) 



X (s'l [e'^" - 1] \s)S{E, +Eg- E,, - Eg,) . (23) 
Using the identity S{x) — l/(27r) J^^dt e"* and employing (sje*^'^*/'' = (sje'^j'''*/'', the above expression becomes 



ft2 



dp f°° dp 



2tt 



— 00 " J —00 



2tt 



vip)v*{p)f,,,{p)f*{p) / dt e^ii^^'^.'^'/'-^ii 



e "B ^ ' — 1 



e "f3 ^ ' — 1 



» (24) 



with ((. ■ ■)) — J2s ^ ^"^"{sl . ■ . \s)/J2s" ^ the average over phonon states. This average can be evaluated for 

q q' employing Glauber's theorem^^. which yields 



- 1 



e "B ^ ' — 1 



-^^^-^'(("(O)')) -T^P^duitf)) 4-PP«"(0)u(t))> 



with 



2p,NsUjQ 



[(1 + neinwQ)) e-'"«* + nB{niOQ)e''^^' 



(25) 
(26) 



the correlation function of the displacement field where pL is the mass of the unit cell of the lattice and Ng 
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is the number of unit cells. 

As in I and II we use for calculational convenience a 
bulk Debye model for the longitudinal acoustic phonon, 
although it is less justified for the high energy part of the 
spectrum which also enters our calculation. Sums over 
phonon momenta are thus replaced by 



3Ns 

, ,3 



duj LO 



(28) 



In terms of the dimensionless variables 

LJ hujJJ 



X 



and T — ujot , 



ujd ksTs 
the phonon correlation function becomes 



(29) 



((u(O)u(t))) 



3h 



da; a; 



1 - e 



— Sx 



- 1 



(30) 



Hence, for the transition probability per unit time we 
obtain 



4a2 foo 1 />oo J ~ foo 
B J-oo ^'^ J-oo ^'^ J-oo 



h-^LUoai .l^r^. 27r 



where 



dxx 



and 7 



3h 



(31) 



(32) 



The transition probability (|5T|) contains two Debye- 
Waller factors, exp(— 7p^q(0)/2) and exp(— 7p^g(0)/2), 
governing the reduction of the transition probability as 
a function of the surface temperature. It also contains 
phonon processes of all orders as can be most easily seen 
from the Taylor expansion 



o7PP'3(t) 



1 + Jppq{T) 



(33) 



Clearly, the second term on the right hand side repre- 
sents the one-phonon and the third term the two-phonon 
process. From I we know that two-phonon processes 
are much weaker than one-phonon processes. We ex- 
pect therefore lower order phonon processes to dominate 
their higher order corrections, so that the expansion p3p 
converges quickly. 

As higher order phonon processes are small compared 
to lower order processes, we take, for a given A^^/ only 
the leading term of exp{jppq{T)) into account, that is, 
the lowest order phonon process that enables a transition 
between the states q and q' . The Fourier transformation 
of powers of (/(r), however, required when p3p is used 
in pil) . cannot be evaluated in closed from, making it 
necessary to construct an approximation for q{T). 

To derive an approximation for (/(r) subject to the con- 
straint 



dre*'^55'^g"(r) = for|A,,.|>n 



(34) 



which states that an n-phonon process yields a non- 
vanishing transition probability only for —n < A < n 
and vanishes otherwise, we split qir) = 9*' (t) -I- 9' (t) into 



I 

a contribution arising from spontaneous phonon emission 
g''(r) and a contribution from induced phonon emission 
or absorption q^ir). They are respectively given by 



q'{T)^ [ dx xe-''''' and q\T) = 2 [ dx 
Jo Jo 



cos(xr) 



1 ■ 

(35) 



(36) 



The former can be evaluated giving 

„ , , cos T — 1 T cos T — sin r sin t 
Q (r) = 5 + * 2 + 

For (/'(t) we need to find an approximation. For that 
purpose we look at the Fourier transform of q^ 

.^Ar„^^_^ _ J ^TT for - 1 < A < 1 



dTe*^"q*(r) 



,5|A 





else 



(37) 



Expanding the Fourier transform in terms of |A| 



2tt- 



;5|A| _ 1 



2tt 



\A\ + -S\A\' + OiS') 



(38) 

yields a high-temperature approximation which con- 
verges quickly for the temperatures we are interested in 
and guarantees at the same time that the one-phonon 
contribution can only bridge energy differences up to 
\Aqq' 1 = 1. Applying the inverse transformation gives 

'2 S\ sin(T) 1 



-1 



6 



6\ cos(t) 6 sin(r) 



0{S'), (39) 
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which satisfies Eq. (p4)) . Using this approximation the 
Fourier transform of powers of (/(r) can be done analyti- 
cally. 

As the n-phonon process gives a vanishing transition 
probability at |A| — n, we take the maximum of the 



n-phonon and the n + 1-phonon process to obtain a bet- 
ter approximation in the vicinity of |A| = n. Then the 
Fourier transformation of exp(7ppq(r)) in leading non- 
vanishing order is given by 



f°° dTe'^'^+'^'P'^'''''^^ w (ma.x{An,An+i) for n - 1 < A < n 

\max(_B„,B„+i) for - n < A < -n -|- 1 ' 



(40) 



where 



- ^ -^^^i: i G) ihh T il ^ r ir^ '-^i^ ^ 



and 



fc=0 ^ ' j=0 



(42) 



Using the approximation given by Eq. (PO)) allows an 
efficient numerical evaluation of the transition probabil- 
ities (IHTI) . 

Equations p9ll42p are first order in 6. For the mate- 
rials and temperatures we are interested in this is suffi- 
cient. Note, however, that the expansion can be contin- 
ued to higher orders in S. The Fourier transformation 
of exp(7ppg(r)) is then still a polynomial in A and thus 
amenable for numerical calculations. 



IV. RESULTS 

We now use the multi-phonon transition probability 
to study the electron kinetics in front of a CaO, AI2O3, 
and Si02 surface. They all have three-phonon deep sur- 
face potentials, that is, the energy difference of the two 
lowest image states is between two and three Debye ener- 
gies. The material parameters required for the numerical 
computation are summarized in Table HI All numerical 
results were obtained for these parameters. Where indi- 
cated, we varied the Debye temperature to simulate dif- 
ferent potential depths. Furthermore, the multi-phonon 
calculation for CaO, AI2O3, and Si02 is compared to the 
one- and two-phonon calculations from I and II which are 
applicable to graphite and MgO. 

A. Desorption 

To judge the quality of the nonlinear multi-phonon ap- 
proximation derived in the previous section, we first com- 
pare in Fig. [T] the inverse desorption time obtained from 



I 

TABLE I. Material parameters for the numerical results. 

CaO AI2O3 Si02 

Debye temperature To 648 K 980 K 470 K 

Dielectric constant £3 12.01 9.9 3.78 

TO-phonon frequency hcuT 41 meV 79 meV 133 meV 

it with the inverse desorption time obtained from our pre- 
vious one- and two-phonon approximation. Shown is the 
dependence of on the Debye temperature Tjy which 
is tuned to vary the potential depth. The dimensionless 
inverse temperature S, as defined in Eq. (I^Hl) . is kept 
constant to keep the level of phonon excitation the same 
while the Debye temperature is varied. 

The nonlinear multi-phonon approximation can be of 
course only compared with the two-phonon approxima- 
tion in the range of Debye temperatures for which the 
potential is two-phonon deep. Calculated in the multi- 
phonon approximation t^^ changes very little over the 
range of two-phonon depth, but shows steep jumps at 
the threshold to one- and three-phonon depth. For 
calculated in the two-phonon approximation, which is 
based on an iteration of the T matrix with the nonlinear 
electron-phonon coupling, these thresholds are washed 
out by the resonances. Nevertheless is on the same 
order of magnitude in both approximations. The main 
effect of the neglected resonances is the rounding off of 
the drops at the thresholds. As the steps, an artefact 
of taking only nonlinear multi-phonon processes into ac- 
count, are less steep for deeper potentials, the nonlinear 
multi-phonon approximation might be even more appro- 
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FIG. 1. (Color online) Inverse desorption time t^^ as a 
function of the Debye temperature Td for 5 — 1 (surface 
temperature Ts — Td/S) calculated in the nonlinear multi- 
phonon approximation (solid red line) and the two-phonon ap- 
proximation from I for two-phonon depth (dashed blue line). 
The surface potential is one-phonon deep for Td > 20007^, 
two-phonon deep for 2000i<r > Td > 1000/C, three-phonon 
deep for lOOOJf > Td > 666A" and four-phonon deep for 
666/C > Td > 500K. Data for Td = 980K apply to AI2O3 
(thin vertical line). 



priate for deeper potentials. The thin dotted vertical line 
in Fig. [1] corresponds to AI2O3. Unfortunately the po- 
tential depth is just below the two-phonon three-phonon 
threshold, so that the value for is most likely under- 
estimated. 

We now move on to the study of the dependence of 
on the surface temperature. Figure [2] shows the in- 
verse desorption time t~-^ as a function of the surface 
temperature for graphite, MgO, AI2O3, CaO, and Si02. 
For graphite and MgO, both two-phonon deep, t~-^ was 
calculated in the two-phonon approximation, for AI2O3, 
CaO and Si02, all of them three-phonon deep, the non- 
linear multi-phonon approximation has been used. For 
all materials t^^ increases significantly with the surface 
temperature. 

Comparing in Fig.[2]T~^ for AI2O3, CaO, and Si02, we 
notice that increases with decreasing Es (see Table |T| 
in accordance with the fact that a smaller implies a less 
deep surface potential and thus a faster desorption. From 
Fig. [5] we also see that for high surface temperatures des- 
orption from the two-phonon deep potentials of graphite 
and MgO is quicker than from the three-phonon deep po- 
tentials of AI2O3, CaO and Si02 as expected. For low 
surface temperature, however, for graphite decreases 
much steeper than for the other materials. This might be 
due to the high Debye temperature of graphite, so that 
for room temperature the dimensionless inverse temper- 
ature S — Td/Ts which controls phonon excitation is 
already in the low temperature regime where downwards 
transitions due to spontaneous phonon emission remain 
constant while upwards transitions are extremely tem- 
perature dependent, causing the desorption time to be 
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FIG. 2. (Color online) Inverse desorption time ^ as a func- 
tion of the surface temperature Ts for graphite, MgO, CaO, 
AI2O3, and SiOz. 



equally temperature dependent. This peculiarity leads 
to the surprising fact that at low temperatures desorp- 
tion from two-phonon deep potentials can be in some 
cases slower than desorption from three-phonon deep po- 
tentials. 

Figure [2] also gives insight into the validity of the high 
temperature expansion (|39[) . Up to first order in b, it is 
valid for (5 < 3. This corresponds to a surface temper- 
ature of 300K for AI2O3 and 225K for CaO. The small 
upwards bends at these temperatures indicate that for 
lower surface temperatures the expansion given by Eq. 
([55)1 should be continued to higher orders in (5. 



B. Sticking 

In II we found that one-phonon processes give much 
higher contributions to the sticking coefhcient than two- 
phonon processes. For this reason we calculate the tran- 
sition probabilities for continuum-bound state transitions 
only in the one-phonon approximation. The effect of 
multi-phonon processes with regard to sticking lies in the 
relaxation from the state the electron is initially trapped 
to the lowest bound state. This is captured by the ki- 
netic sticking coefficient. Before we address this question 
in more detail we take a look at the prompt sticking co- 
efficient. 

The prompt sticking coefhcient for AI2O3, CaO and 
Si02 is presented in Fig. |31 First we consider the prompt 
energy-resolved sticking coefficient shown in the inset. 
Note that the quadratic phonon dispersion of the Debye 
model translates into an energy-resolved sticking coef- 
ficient which, apart from the discontinuities, is propor- 
tional to the electron energy. The steep jumps in the 
energy-resolved sticking coefficient reflect level accessibil- 
ity. When the energy difference between the approach- 
ing electron and a bound state exceeds the Debye en- 
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FIG. 3. (Color online) Prompt energy- averaged sticking co- 
efficient for CaO, AI2O3, and Si02 as a function of the mean 
energy of the electron at a surface temperature of Ts — 300K. 
Inset: Prompt energy-resolved sticking coefficient for AI2O3 
and Si02 as a function of the electron energy for = SOO-K". 
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FIG. 4. (Color online) Prompt (full line) and kinetic (dashed 
line) energy averaged sticking coefficient for Si02 as a function 
of the mean energy of the electron and the surface tempera- 
ture Ts. 



ergy, one-phonon processes no longer enable sticking to 
this level. Since the lowest two bound states of the im- 
age potential of AI2O3 and Si02 have energies e„ < — 1, 
they cannot be reached by onc-phonon processes from the 
continuum. Thus, the lowest bound state contributing to 
prompt sticking is the third bound state. Due to the dif- 
ferences in the Debye energy, the energy-resolved sticking 
coefficient for Si02 is larger and increases faster for low 
electron energies than the sticking coefficients for AI2O3 
and CaO (not explicitly shown). Compared to AI2O3 
and CaO the energy-resolved sticking coefficient for Si02 
is thus strongly peaked at low electron energies. As a 
result, the energy averaged-sticking coefficient shown in 
the main panel of Fig. [3] is much larger for Si02 than for 
AI2O3 and CaO. 

Figure 0] shows the prompt and kinetic energy aver- 
aged sticking coefficient for Si02 as a function of the 
mean electron energy and the surface temperature. The 
prompt sticking coefficient increases slightly with tem- 
perature due to the increased contribution of induced 
phonon emission responsible for continuum-bound state 
transitions. The kinetic sticking coefficient is smaller 
than the prompt sticking coefficient by four to five or- 
ders of magnitude and decreases with temperature as a 
higher surface temperature favors quick transitions back 
into the continuum after initial trapping. 

Depending on whether transitions from the upper 
bound states to the lowest state or to the continuum 
are more hkely, the electron trickles through after initial 
trapping or desorbs before relaxing to the lowest bound 
state. For the three-phonon deep surface potentials of 
AI2O3, CaO and Si02 trickling through is suppressed 
leading to a considerable reduction of the kinetic com- 
pared to the prompt sticking coefficient. 



C. Electron kinetics 

So far we have calculated from the kinetic rate equa- 
tion ^ the prompt and kinetic sticking coefficients and 
the desorption time. The rate equation contains how- 
ever more information. For a specified electron influx or 
initial condition, the time-evolution of the bound state 
occupancy and the energy resolution of the desorption 
flux can be calculated as well. 

To address the first question, we plot in Fig.[5]the time 
evolution of the bound state occupancy. Our aim is to 
identify the stages of physisorption and to relate them 
to the eigenvalues of the matrix Tnm- The situation we 
are considering is a three-phonon deep potential, that 
is, —3 < A12 < —2, with the second bound state lying 
more than one Debye energy below the continuum, that 
is, £2 < —1, so that electron trapping, due to one-phonon 
processes, can occur only in the third and higher bound 
states. The bound states i > 2, which we call upper 
bound states, are linked by one-phonon transitions. 

To obtain the time evolution of the bound state oc- 
cupancy after trapping of an electron at t = 0, we solve 
the rate equation with the initial condition for the bound 
state occupancy, ni{0) = TtWik, for a specified k, which 
is the probability that the electron is trapped in state 
i. In Fig. [SJ k corresponds to an electron energy of 
E = 0.05eV. Due to the high electron energy the third 
state cannot be reached by a one-phonon process. Thus 
trapping occurs in the fourth and higher bound states. 

In a first stage after trapping of the electron, the fast 
one-phonon transitions between the upper bound states 
dominate the electron kinetics. Due to trapping of the 
electron in the fourth and higher bound states, the oc- 
cupancy of the upper bound states is out of equilibrium. 
Over the time scale set by A^^, the inverse of the third 
eigenvalue of T, the occupancy in the upper bound states 
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FIG. 5. (Color online) Time evolution of the bound state 
occupancy of a single electron trapped at f = in the upper 
bound states of AI2O3 at Ts — 300K. The thin vertical lines 
correspond to A^^, Aj~^, and A^^, respectively, where — Ai are 
the three lowest eigenvalues of the matrix T. 
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FIG. 6. (Color online) Energy-resolved desorption flux at 
t — 10~^s for an electron tapped at t = in the upper bound 
states under the same conditions as in Fig.[5l For AI2O3 one- 
phonon transitions between bound and continuum states are 
only possible from the third and higher bound state. The 
small numbers give the bound state from which transitions to 
the continuum are no longer possible at the respective energy. 



relaxes towards its equilibrium value. The electron trick- 
les through from the fourth and higher bound states to 
the second bound state, as can be seen from the increase 
in the occupancy of the second bound state 712 and the 
reduction of the occupancy of the third and higher bound 
states ng>3. 

Then the strong one-phonon transitions between the 
upper bound states and the continuum, occurring over 
the time scale set by Aj~^, empty the upper bound states. 
The weak multi-phonon transitions from the upper states 
to the lowest bound state are only a small perturbation 
to the electron kinetics in the upper bound states so that 
Aj~^ corresponds to the desorption time for the system of 
the upper bound states without the lowest bound state. 

Until the upper bound states are emptied, a small frac- 
tion of the occupancy reaches the lowest bound state as 
can be seen from the discrepancy of the initial occupancy 
of the upper bound states nq>2 and the maximum occu- 
pancy of ni . This difference corresponds to the reduction 
of the kinetic with respect to the prompt sticking coeffi- 
cient in Fig. 21 The lowest bound state remains occupied 
for a much longer time, until desorption takes place at 
times on the order of Te = . 

Figure |6] finally shows the energy-resolved desorption 
flux a.t t = 10~^s (given by Eq. ^) after trapping 
of the electron under the same conditions as in Fig. [S] 
The final transition that sets the electron free is a one- 
phonon transition from one of the upper bound states 
to the continuum. From each bound state i, one-phonon 
transitions are only possible to continuum states with 
an energy E < Ei + hwo- Hence, the energy-resolved 
desorption flux exhibits the same discontinuities as the 
energy-resolved sticking coefficient shown Fig. [31 located 
at electron energies for which one-phonon transitions be- 
tween bound states and the continuum cease to be oper- 



ational. 



V. TWO-STATE SYSTEM: DISCUSSION 

To clarify the generic behavior of electron physisorp- 
tion at dielectric surfaces, and to put the results pre- 
sented in the previous section and in I and II into per- 
spective, we study a simple model of two bound states 
coupled to a continuum of states. Electron physisorp- 
tion occurs in the image potential which supports a deep 
lowest bound state, well separated from a group of rela- 
tively closely packed upper bound states. Since the upper 
bound states are strongly coupled by one-phonon pro- 
cesses they can be subsumed under an effective upper 
bound state. The effective state is then weakly coupled 
to the lowest bound state via multi-phonon processes and 
strongly coupled to the continuum via one-phonon pro- 
cesses. 

The left panel of Fig. |7l schematically shows the system 
of the two surface states. Gearing towards deep poten- 
tials, we include only transitions between the two bound 
states and between the upper state and the continuum. 
The matrix T defined by Eq. (JT]) reads for this system 



W21 -Wi2-Wc2, 



(43) 



where W12 and W21 are the transition probabilities from 
the second to the first bound state and vice versa, and 
Wc2 is the transition probability from the second bound 
state to the continuum. For the two-state model the 
eigenvalues — and the right and left eigenvectors e^'^-' 
and e'^'^) can be calculated analytically. The eigenvalues 
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FIG. 7. (Color online) Left panel: Schematic drawing of the 
two-state model discussed in the main text. Middle panel: 
Physisorption scenario of type A. A trapped electron has a 
high chance to drop to the bottom, then it revolves between 
the two bound states, until it desorbs. Right panel: Ph- 
ysisorption scenario of type B. Due to a relaxation bottleneck 
the electron is unlikely to drop to the lowest state. Transi- 
tions that the electron makes once per temporary trapping 
event are represented by a thin line, a bold line represents 
transitions made more than once, and dashed lines represent 
transitions that are made with a very low probability. 



are given by 



-Ao,i = - i {W21 + W12 + W,2) 



± 2 V (^21 + W12 + Wc2) - m^2lWc2 ■ (44) 

Parameters of physical interest can be also obtained 
analytically. The desorption time, for instance, is the 
negative inverse of the lowest eigenvalue, Te = ■ Since 
only the upper bound state can be reached from the con- 
tinuum, prompt sticking arises solely from trapping in 
the upper bound state. The prompt sticking coefficient 
is thus given by s^™™^* = TtW2k- In the two-state model, 
the kinetic sticking coefficient is moreover related to the 
prompt sticking coefficient by s^'" = 62"'' sP™™^*. Hence, 
the probability for the electron to trickle through from 
the upper to the lower bound state is 62^^^ 

For many dielectrics the weakest transitions are from 
the lowest bound state to the upper bound states. They 
are typically triggered by more than two phonons. To 
mimic this situation within the two-state model we set 
W21 ^ Wi2,Wc2- The inverse of the desorption time 
becomes in this limit 



W, 



c2 



W12 -f Wc 



-W2I 



(45) 



and the ratio between kinetic and prompt sticking coef- 
ficient becomes 



„kin 



prompt 



W12 + Wc2 



(46) 



The physical behavior of the two-state model depends 
therefore on the ratio between Wi2 and Wc2 and thus 
on the potential depth and the surface temperature. 
Two extreme cases are possible and represent differ- 
ent physisorption scenarios. For Wc2 ^ W^i2, t~'^ = 



{W2i/Wi2)Wc2, which, using detailed balance, can be 
brought into the Arrhenius form = e~^^^^~^^^Wc2- 
Kinetic and prompt sticking coefficients coincide more- 
over in this parameter range. Hence, an electron trapped 
in the upper state drops to the lowest state before des- 
orption. Desorption from the lowest state occurs then via 
a cascade, that is, a series of fast transitions 1 — > 2 
1 — 2 — > 1 . . . until eventually the transition 2 — cont. 
removes the electron from the bound states. The just 
described physisorption scenario, which we call type A 
scenario, is illustrated in the middle panel of [T] Recall- 
ing that the upper level stands for a manifold of strongly 
coupled bound states, it resembles the physisorption of 
neutral particles via cascades originally proposed and in- 
vestigated by Gortel and coworkers.— 

In the other limit, W12 ^ Wc2- The inverse of the 
desorption time and the ratio between prompt and ki- 
netic sticking coefficient are then given by = W21 
and gkin/gprompt ^ Wi2/Wc2, respectively. The ph- 
ysisorption scenario is now dramatically different from 
the one discussed before because the desorption time is 
solely determined by the transition probability from the 
lower to the upper bound state. As a result, desorption 
does not occur via a cascade, but as a one-way process 
1 ^ 2 — !> cont., where the second transition is so fast that 
it basically does not affect the desorption time. Hence, 
in this scenario, which we call type B, the upper bound 
state can be considered as de-facto belonging to the con- 
tinuum and desorption as basically equivalent to desorp- 
tion from a single deep state. For sticking, the type B 
scenario exhibits moreover a relaxation bottleneck. An 
electron trapped in the upper state is very unlikely to 
drop to the lowest bound state, as schematically shown 
in the right panel of Fig. [T] 

Within the limits set by the model for the electron- 
surface interaction introduced in I and briefly recalled 
in section IhH the two-state model contains the essential 
physics of electron physisorption. For potentials with 
£2 > —1, where a direct one-phonon transition from 
the second bound state to the continuum is possible the 
two-state model can be applied directly. Calculating the 
desorption time within the two-state model shows very 
good agreement with the results for graphite obtained 
in I. For the ratio between kinetic and prompt sticking 
coefficient skmy^prompt^ which is given in the two-state 

model by e2°\ the agreement is less good but qualita- 
tively correct, reproducing, for instance, the tempera- 
ture dependent transition between type A and type B. 
For potentials with £2 < — 1 no one-phonon process from 
the second bound state to the continuum is possible and 
the two-state model cannot be applied directly. For ph- 
ysisorption of type B, however, the electron kinetics in 
the upper bound states is only marginally perturbed by 
transitions to and from the lowest bound state. The time 
it takes an electron to get from the second bound state 
into the continuum is then the desorption time of the 
system of the upper bound states alone, that is, the neg- 
ative inverse of the smallest eigenvalue — Aq^ of the ma- 
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trix T"Pj which is the matrix Tnm defined in Eq. ^ 
with n,m > 1. In the two-state model Aq^ can be re- 
garded as an effective transition rate between the second 
state and the continuum. Hence, to apply the two-state 
model with potentials where the second bound state does 
not couple by one-phonon processes to the continuum we 
simply replace in Eqs. P5|) and (^5)) Wc2 by Aq^. 

Let us finally look at the results obtained in the pre- 
vious section and in I and II from the perspective of the 
two-state model. In a one-phonon deep potential, the 
transitions from the upper bound states to the lowest 
bound state and from the upper bound states to the con- 
tinuum are enabled by one-phonon processes. In this case 
the downward transitions are always more likely than the 
upward transitions so that one-phonon deep potentials 
give always rise to physisorption of type A. Hence, they 
show no relaxation bottleneck and prompt and kinetic 
sticking coefficient coincide. Two- or more-phonon deep 
potentials can either lead to physisorption of type A or 
type B, depending on the surface temperature. In this 
case one-phonon transitions from the upper bound states 
to the continuum compete with multi-phonon transitions 
from the upper bound states to the lowest bound state. 
As a transition from the upper states to the continuum 
requires phonon absorption, proportional to ub, while 
a transition from the upper state to the lower requires 
phonon emission, proportional to I + ub, we expect that 
for sufficiently low temperature physisorption is always of 
type A, even for multi-phonon deep potentials. For suffi- 
ciently high temperatures all two- or more-phonon deep 
potentials are however of type B. In this case a relaxation 
bottleneck results in the discrepancy between prompt 
and kinetic sticking coefficient (see Fig. 2]). The elec- 
tron kinetics is primarily determined by the one-phonon 
transitions among the upper states (see FigO. The tem- 
perature at which type A merges into type B depends 
on the potential depth and the Debye temperature. For 
room temperature the three-phonon deep potentials of 
AI2O3, CaO and Si02 and the two-phonon deep poten- 
tial of MgO are all of type B. The crossover between type 
A and type B occurs for the two-phonon deep potential of 
graphite at room temperature (see Fig. 5 of our previous 
work II). 



VI. CONCLUSIONS 

Within a simplified one-dimensional model for the 
polarization-induced interaction between an external 
electron and a dielectric surface with a sufficiently large 
energy gap and a sufficiently negative electron affinity, 
we investigated phonon-induced adsorption and desorp- 
tion of an electron at a dielectric surface. The required 
electron energy relaxation, inducing transitions between 
the eigenstates of the surface potential, which we approx- 
imated by a recoil-corrected image potential, is due to the 



coupling to an acoustic bulk phonon. 

The majority of dielectrics of interest have a surface 
potential that is three- or more-phonon deep, that is, the 
energy difference between the two lowest bound states is 
more than two Debye energies of the bulk phonon. In 
our previous work j^^i^^ we took multi-phonon processes 
into account using a T matrix approach, which is however 
only feasible for one- and two-phonon deep potentials, as 
it is, for instance, the case for graphite. To overcome this 
limitation we derived in this work a non-perturbative ex- 
pression for the multi-phonon transition probability aris- 
ing solely from the nonlinearity of the electron-phonon 
interaction. In view of our previous results for one- and 
two-phonon deep potentials we expect this approxima- 
tion to give an acceptable correct order of magnitude 
estimate for the multi-phonon transition probability in- 
volving more than two phonons, despite the neglect of 
resonant processes stemming from the iteration of the T 
matrix. 

We presented numerical results for the electron des- 
orption time for graphite, MgO, CaO, AI2O3, and Si02 
and the prompt and kinetic energy-resolved and energy- 
averaged electron sticking coefficient for CaO, AI2O3, 
and Si02. In addition, we calculated the energy- resolved 
desorption flux and investigated the time evolution of 
the bound state occupancy after initial trapping of an 
electron, revealing the characteristic stages of electron 
physisorption: initial trapping, relaxation in the upper 
bound states, trickling through to the lowest bound state, 
and desorption. Ultrafast electron spectroscopy at sur- 
faces with stable image states^"— should be able to resolve 
these different stages experimentally. 

Using a simple two-state model, we finally identified 
two vastly different scenarios of electron physisorption, 
depending on potential depth and surface temperature, 
and put our results, including the ones of our previous 
work^'^'^'*, into perspective. For almost all dielectrics of 
practical interest, the trapped electron has only for very 
low temperatures, well below room temperature, a signif- 
icant chance to trickle through to the lowest bound state. 
The desorption process in this case would then proceed 
via a cascade between the first and second bound state, 
until it eventually makes a transition from there to the 
continuum. The shallow bound states, albeit important 
for adsorption, play a minor role for desorption. The 
second bound state is the most important one. It is a 
relay state. At room temperature, however, a relaxation 
bottleneck prevents the trapped electron from falling to 
the lowest bound state. The electron physisorption ki- 
netics is thus dominated by fast one-phonon transitions 
in the upper bound states. Only a small fraction of the 
electron trickles through to the lowest bound state and 
resides there for a very long time until it makes a one-step 
desorbing transition to the continuum. 
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